Plasmon-phonon coupling in graphene 
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Collective excitations of coupled electron-phonon systems are calculated for both monolayer and 
bilayer graphene, taking into account the non-pert urbative Coulomb coupling between electronic 
excitations in graphene and the substrate longitudinal optical phonon modes. We find that the 
plasmon-phonon coupling in monolayer graphene is strong at all densities, but in bilayer graphene 
the coupling is significant only at high densities satisfying the resonant condition cj p i « uj p u . The dif- 
ference arises from the peculiar screening properties associated with chirality of graphene. Plasmon- 
phonon coupling explains the measured quasi- linear plasmon dispersion in the long wavelength limit, 
thus resolving a puzzle in the experimental observations. 
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A plasmon is a collective mode of charge-density os- 
cillation in the free-carrier system, which is present 
both in classical and quantum plasmas. Studying the 
collective plasmon excitation in the electron gas has 
been among the very first theoretical quantum mechani- 
cal many-body problems studied in solid-state physics. 
The collective plasmon modes of monolayer graphene 
(MLG) have been extensively studied theoretically^— 
and experimentally^ - — . Recent discovery of bilayer 
graphene (BLG) has also led to a number of theoreti- 
cal descriptions of plasmon modes in BL G 12 i 13 . Even 
though the long-wavelength plasmon frequency of MLG 
is explicitly nonclassical (i.e., the plasmon frequency is 
necessarily quantum with H appearing manifestly in the 
long- wavelength plasma frequency 3 -) , its wave vector dis- 
persion is given by classical electrodynamics, i.e., uo p (q) oc 
y/q. Note that the quadratic band-dispersion of BLG 
makes the leading order long- wavelength plasmon disper- 
sion explicitly classical with the same ~ ^fq dependence. 
However, at finite g, away from the long- wavelength limit, 
there are several corrections to the plasmon dispersion 
(j(q) arising from nonlocal finite wave vector response, 
finite-temperature thermal corrections, many-body ef- 
fects, local-field corrections, and other mechanisms rel- 
evant to the specific electron system. 

Since the plasmon dispersion relation is exactly known 
at long wavelengths (q — »• 0) where the f-sum rule 
arising from particle conservation fixes the plasma fre- 
quency, it is surprising that the measured graphene plas- 
mon dispersion in the long wavelength limit deviates 
from the classical dispersion (u q ~ y/q) and shows a 
rather linear dispersion^^. In a recent experiment 9 , the 
strongly coupled plasmon-phonon mode dispersion has 
been measured by the angle-resolved reflection electron- 
energy-loss spectroscopy and it is found that the dis- 
crepancy arises from electron-phonon coupling. In epi- 
taxial graphene the substrate (i.e. SiC) is a highly 
polar material. In general, carriers in polar materi- 
als couple with the longitudinal optical (LO) phonons 
of the system via the long-range Frohlich interaction. 
However, the surface optical (SO) phonon is a well- 
characterized surface property of polar semiconductors, 



and it is possible that carriers in graphene layer cou- 
ple to the SO-phonons of the underlying substrate lat- 
tice via the long-range polar Frohlich coupling 14,15 . For 
isotropic media the frequency of SO phonons uoso is re- 
lated to the transverse optical (TO) bulk phonon ujto as 

uso/uto = V (eo + l)/(e<x> + ljr 4 -"— where e (eoo) is 
the static (high frequency) dielectric constant. Note that 
the bulk longitudinal optical phonons ujlo and ujto are 
connected with the dielectric consta nts by the Lyddane- 
Sachs- Teller relation uolo/^to = \/ e o/ e °o- 

The electron-phonon coupling is the macroscopic cou- 
pling of the electronic collective modes (plasmons) to 
the optical phonons. The mode coupling phenomenon, 
which hybridizes the collective plasmon modes of the 
electron gas with the optical-phonon modes of the lat- 
tice, gives rise to the coupled plasmon-phonon modes (the 
hybrid modes) which have been extensively studied^— 
both experimentally and theoretically in bulk and 2D 
electron systems. The electron-phonon interaction leads 
to many-body renormalization of the single-particle free 
carrier properties 21,22 and also affects the transport 
propertie s 23 ! 24 . A good understanding of electron- 
phonon coupling is thus important in developing quanti- 
tative theories for many different experimental studies in 
graphene. 

In this paper we calculate the coupled plasmon-SO 
phonon modes of epitaxial graphene (or graphene on a 
polar substrate such as Si02, SiC, or Hf02). Our most 
significant finding is that in MLG plasmon-phonon mode 
coupling effect is strong at all electronic densities due 
to the singular behavior in the screening function aris- 
ing from chirality 1 . By contrast, for BLG, the plasmon- 
phonon coupling is significant only at high carrier den- 
sities. We also find that at low densities, when the cou- 
pling is weak and the coupled phonon-like mode (gapped 
mode) lies in the interband electron-hole continuum, the 
energy of phonon-like mode decreases in the long wave- 
length limit due to the coupling of the phonon mode to 
the interband single particle excitation, which arises from 
the enhanced BLG backscattering 25 . However, at high 
densities, when the plasmon-phonon mode coupling is 
strong, the phonon-like mode frequency increases linearly 
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FIG. 1. Effective dynamical interaction V e ff calculated in 
RPA. Dashed (wiggly) lines represent the SO-phonon medi- 
ated (Coulomb) electron electron interaction V p h (v c ), and 
the bubble the irreducible polarizability U(q,cj). 



with wavevector, as in MLG. 

We first present our model for plasmon-phonon cou- 
pling, which consists of a two-dimensional electron gas 
coupled to dispersionless SO-phonons at zero tempera- 
ture. For MLG, we have a system of Dirac fermions with 
linear dispersion, while BLG have a parabolic dispersion 
around the Dirac point. Due to the presence of the long 
range electron-phonon coupling, electrons interact among 
themselves through the Coulomb interaction and through 
virtual- SO-phonon exchange via the Frohlich interaction. 
The electron- SO phonon interaction is given by 
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kq ss' 



J 
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Cks (b q + bl^j , (1) 



where c\ s is the electron (s = +1) or hole (s = —1) 
creation operator, b q and b q are creation and destruction 
operators of surface phonon, and the interaction matrix 
element M| s is defined by 



Mg =M (q)Fi +q F l 



s'ki 



(2) 



where F sk is the chiral spinor and given by F^ k = 
(s,e~ iek )/V2 with s = ±1, k = tarn- 1 (k y / k x ) for 
MLG^ and k = 2 tan" 1 (k y / 'k x ) for BLG 2 ^. We also 
have 



Q ^ 



eo 



(3) 



where d is the separation distance between graphene layer 
and substrate. Then the matrix element of SO-phonon 
mediated electron-electron interaction is dependent on 
both wave vector and frequency and give by, 



v ph {q,u) = [M (g)] 2 A)M, 



(4) 



where the unperturbed SO-phonon propagator is given 
by D (uj) = 2uj S o/(oj 2 - oj 2 so ). 

The total effective electron-electron interaction is ob- 
tained in RPA 19 by summing all the bare bubble dia- 
grams (see Fig. [I]), 



= v c (q) 

1 - [v c (q) +v ph (q,uj)}no(q,uj) e t (q,u) 



v c (q) +v ph (q,u) 



(5) 

where v c (q) = 2ne 2 /e 00 q is the electron-electron 
Coulomb interaction and Ho(q : oj) is the complex irre- 
ducible polarizability of either the monolayer— or the bi- 
layer system^ 3 given by the bare bubble diagram. The 
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FIG. 2. Calculated plasmon-phonon coupled mode uj± disper- 
sions in MLG as a function of the wave vector q for two differ- 
ent densities (a) n — 10 13 cm~ 2 and (b) n — 10 12 cm -2 . The 
plasmon dispersion (uj q ) without the electron-phonon cou- 
pling is shown by the dashed line. Two dotted horizontal 
lines represent the frequencies of the uncoupled SO (top) and 
TO (bottom) phonon modes, respectively. 



total dielectric function within RPA contains contribu- 
tions both from electrons and SO-phonons: 



e t {q,u>) = 1 U {q,uj)- 



where 



ae~ 2qd 
l-ae- 2 « d -uj 2 /uj 2 so J (6) 



(7) 



1 e + l 



The collective mode dispersion is given by the zeros of 
the complex total dielectric function: e t (q,uj) — 0. 

Let us first focus on the collective modes of MLG. In 
the long wavelength limit (q —> 0) we get the following 
coupled uj± collective modes 



u+(q) uj S o 
u-(q) = (l- 



-2qd . 



uso 
ae~ 2qd ) uo n 



(8a) 
(8b) 



where uo 2 = 2e 2 EFq/e OQ (Ep = Fermi energy) is the 
plasmon mode dispersion of an uncoupled system in the 
long wave length limit. As q — >• the phonon-like mode 
uu+ is located above uoso and increases linearly, and the 
plasmon-like uo- is slightly less than the corresponding 
uncoupled monolayer graphene plasmon mode, uo q . 

In Fig. [2] we show the calculated coupled plasmon- 
phonon collective modes in MLG for two different densi- 
ties. The following parameters are used throughout this 
paper 29 : ujto = 95.0 meV, uso = 116.7 meV, = 6.4, 
eo = 10.0, and d = 5A. As shown in Fig. [2] the mode 
coupling in MLG is strong for all electron densities. In 
ordinary 2D systems or 3D systems the plasmon-phonon 
mode coupling is only significant at densities satisfying 
the resonant condition u q « uoso- However, in MLG the 
plasmon mode exists for all wave vectors due to the singu- 
lar behavior in the polarizability, which leads to strong 
plasmon-phonon coupling. Since the singular behavior 
of the polarizability arise from the suppression of the 
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FIG. 3. The density plots of energy loss function 
(— Im[l/e(g, uj)]) of MLG in (q, uj) space for two different den- 
sities (a) n = 10 13 cra~ 2 and (b) n = 10 12 cra~ 2 . Note that 
the uj-\- mode of the high density result carries most spectral 
weight and its dispersion is almost linear, which is observed in 
the recent experiment 8 . We use a phenomenological damping 
of 0.1 i^F in these results. 



FIG. 4. Calculated plasmon-phonon coupled mode uj± dis- 
persions in bilayer graphene as a function of the wave vec- 
tor q for two different densities (a) n = 10 13 cra -2 and (b) 
n = 10 12 era -2 . The plasmon dispersion (cj g ) without the 
electron-phonon coupling is shown by the dashed line. The 
dotted horizontal line represent the frequency of the uncou- 
pled SO phonon mode. 



back scattering due to the chirality of MLG the strong 
plasmon-phonon coupling is a direct consequence of its 
unique chiral property of MLG. Note that the plasmon- 
like mode uo- in Fig. [2] vanishes at a finite critical wave 
vector, q c c± ujso(^ — ol)/vf, and for q > q c we find only 
the phonon- like mode (gj+) which approaches u q for large 

Q ' 

The dynamical structure factor, S(q,u), which gives 
the spectral weight of the collective modes, is propor- 
tional to the imaginary part of the inverse dielectric func- 
tion (loss function) and given by 



S(q,u>) 



1 



nov c (q) 



Im 



(9) 



For a true collective mode with zero Landau damping 
both Im[e t (q,uj)] and Re[e t (q,uj)] vanish, and the inverse 
dielectric function becomes a delta function with weight 



£Ree t (<2,cj)| 



(10) 



where uji (q) is the collective mode frequency at wave vec- 
tor q. In the long wavelength limit the weight of plasmon- 
like mode can be calculated as 



^(i-«) 3/ V 



W(q)\ u . 

and the weight of phonon-like mode as 



W(q)\ u+ = irawso/2- 



(11) 



(12) 



The spectral weight of uj- mode vanishes as yjq in the 
long wavelength limit, but the weight of cj + mode is finite. 
Thus in the long wavelength limit all spectral weight is 
carried by the phonon-like mode. In Fig.[3]the calculated 
loss function — Im[l/e(g, uj)] is shown in (q,uj) space for 
two different densities (a) n = 10 13 cm~ 2 and (b) n = 
10 12 cm -2 . In the long wavelength limit the phonon-like 
mode has most of the weight. In the intermediate wave 



vector range, however, the plasmon-like mode becomes 
stronger. The weight of the uj- mode vanishes again 
when the plasmon-like mode merges with the electron- 
hole continuum at a critical wave vector and uj- mode 
becomes overdamped by Landau damping. 

Let us now turn our attention to BLG. Just like MLG, 
we again get two hybridized plasmon-phonon modes, one 
(uj-) having a ~ ^fq dispersion and the other (cj+) ex- 
hibiting a gap equal to the SO phonon frequency (ujso) 
in the long wavelength limit. The cj_ mode has the same 
dispersion as in MLG, cu-(q) = (l-ae- 2qd )uj q , which lies 
in the gap between the intraband and interband continua 
and has a spectral weight which goes as ~ y/q. Thus, in 
the long wavelength limit, all the oscillator strength lies 
in the gapped mode cj + . 

However, there are two main differences from MLG, 
i.e., the quadratic energy dispersion and the enhanced 
backscattering due to chirality in BLG 25 , which lead to 
non-trivial differences in the collective mode spectrum. 
These two effects lead to very different behaviour in the 
low and high density limits. To illustrate these effects, 
we plot the collective mode spectrum of bilayer graphene 
at two different densities, (a) n = 10 13 cm -2 (high den- 
sity, Ep > ujso) and (b) n = 10 12 cm -2 (low density, 
Ep < oo so) in Fig.Hl Here, u q is the uncoupled plasmon 
frequency and the shaded regions represent the intraband 
and interband particle-hole continuum. The correspond- 
ing loss functions are plotted in Fig. [5) 

In the high density limit where w q ~ ojso, there is 
a strong plasmon-phonon coupling as evidenced by the 
deviations of cj + from ujso and of uo- from uj qi which gives 
rise to the gapped mode cj + having a linear dispersion 
with a positive slope in the low q limit. At larger q values, 
it approaches the uncoupled plasmon dispersion, as seen 
in Fig. H](a). The uj- mode merges into the continuum 
at a critical wavevector q c , which is much smaller than 
that of the uncoupled mode indicating strong electron- 
phonon coupling. Furthermore, as seen from Fig. [5fa), 
the phonon-like mode cj + carries a much larger spectral 
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FIG. 5. The density plots of energy loss function 
(— Im[l/e(g, cj)]) of bilayer graphene in (q,u) space for two 
different densities (a) n = 10 13 cm~ 2 and (b) n — 10 12 cm~ 2 . 



weight than uj-. 

In the low density limit where u q < uuso, the plasmon- 
phonon coupling is weak and the gapped mode cj + is 
barely affected by the coupling. In addition, the mode 
energy decreases linearly in the long wave length limit, 
as seen in Fig. H(b). The small negative slope is a con- 
sequence of the coupling of the phonon mode to the in- 
terband single particle excitation^ arising from the en- 
hanced backscattering in the system and is a distinct dif- 
ference between the MLG and BLG. Note that when the 
SO phonon mode is pushed into the interband electron- 



hole continuum, the coupled uj+ mode is always Landau- 
damped due to the presence of the interband continuum 
and carries little spectral weight beyond a very small 
range of low q values. The deviation of the plasmon-like 
mode uj- from the uncoupled dispersion is much smaller 
than in the high density limit, further showing that the 
plasmon-phonon coupling is weak in this limit. From 
Fig. [5fb), we find that beyond a small range of low q 
values, the plasmon-like mode carries much more spec- 
tral weight than the phonon-like mode and hence, at low 
densities, the plasmon mode should be easier to detect in 
BLG. 

In summary, we have calculated the dispersion and 
the spectral weight of the coupled plasmon-phonon mode 
of 2D graphene. We find that the mode-coupling effect 
is strong in monolayer graphene at all densities in con- 
trast to the corresponding bilayer graphene, where the 
coupling is only significant at high densities. Since the 
carriers in graphene are strongly coupled to the surface 
optical phone of a polar substrate it is important to un- 
derstand the many-body renormalization of the single- 
particle properties and the transport properties in the 
presence of electron-SO phonon coupling. 
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